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SUMMARY 


(^e  buckling  strengths  of  eleven  axially  compressed  filament- 
wound”* cylinders  with  diameter  to  thickness  ratios  ranging  from 
167  to  643  are  determined  experimentally  and  analytically.  The 
analytical  predictions  of  buckling  loads  are  made  by  using  linear 
anisotropic  shell  theory.  The  results  of  the  compression  tests 
generally  indicate  that  the  cylinders  buckled  at  65  to  857o  of  the 
loads  predicted  by  linear  analysis.  The  composite  moduli  for 
each. of  the  cylinders  are  determined  experimentally  and  are  in 
reasonable  agreement  with  the  theory  of  S.  Tsai.f 


SYMBOLS 


In-plane  stiffness  matrix 


Element  of  characteristic  determinant 


Modal  amplitudes 


^xx’ 


Coupling  stiffness  matrix 

Cos  cpj  also  used  as  a  superscript  or  subscript 
to  indicate  circumferential  layer 

S.  Tsai’s  contiguity  factor 

Modulus  matrix 

Composite  modulus  perpendicular  to  the  fibers 

Composite  modulus  representing  Poisson’s  ratio 
effect 

Composite  modulus  in  the  fiber  direction 

Composite  shear  modulus 

Mean  cylinder  diameter 
Bending  stiffness  matrix 

Strains 

Young’s  modulus  of  filament 
Young’s  modulus  of  matrix 


Composite  Young’s  modulus  perpendicular  to  the 
fibers  =  (1  -  •  ^2l) 


Axial,  circumferential,  and  radial  displacements 
Axial,  angular,  and  radial  coordinates 

Cartesian  coordinate  system 
o 

j 

In -plane  strains 

Specific  gravity  of  filament 

Specific  gravity  of  matrix 

Axial  wave  number 
Poisson's  ratio  of  filament 

Poisson's  ratio  of  matrix 

Major  Poisson's  ratio  of  composite 

Minor  Poisson's  ratio  of  composite 


Angle  of  polar  winding 


I.  INTRODUCTION 


Experimental  data  on  buckling  of  f ilament -wound  cylinders 

under  axial  compression  have  been  obtained  by  Young,  Card  and 

Peterson,^  and  Sargent  (see  Ravenhall^) ,  The  cylinders  tested 
had  diameter  to  thickness  ratios  ranging  from  25  to  304.  In 
general,  a  reasonable  comparison  was  achieved  between  the  buck¬ 
ling  loads  predicted  by  small  deflection  orthotropic  shell 
theory  and  the  test  data.  In  addition,  Card  and  Peterson  noted 
the  coupling  between  shear  failure  and  buckling  when  resin 
stresses  are  high  enough  at  the  buckling  load  to  cause  plastic 
flow  in  the  resin. 


This  report  describes  an  experimental  and  analytical  study 
of  the  buckling  strength  of  eleven  filament -wound  cylinders  under 
axial  compression.  The  cylinders  fabricated  for  this  study  have 
diameter  to  thickness  ratios  of  167  to  643  and  consist  of  three 
composite  layers.  Fabrication  procedures  for  the  cylinders  are 
given  in  Appendix  A.  One  layer  is  a  polar  wrap  consisting  of  two 
half-layers  at  angles  of  -Kp  and  -cp  with  respect  to  the  longitudinal 
axis  of  the  cylinder  (0  <  cp  <  9  deg).  The  polar  wrap  is  enclosed 
by  an  inner  and  an  outer  circumferential  wrap.  The  three  layers 
have  the  same  nominal  thickness.  This  combination  is  chosen  be¬ 
cause  it  is  efficient  for  internal  pressure  and  minimizes  the  ef¬ 
fects  of  initial  imperfections  on  buckling  strength  by  providing 

. -  -  .  4,5,  and  6  ... 

circumferential  in-plane  and  bending  stiffening.  Ail 

of  the  cylinders  studied  have  in-plane  and  bending  stiffness 

matrices  of  orthotropic  form.  Also,  the  cylinders  whose  middle 

layers  are  not  oriented  axially  produce  an  anisotropic  coupling 

between  stretching  and  shear  of  the  form  described  by  Reissner 

and  Stavsky.^ 


Analytical  predictions  for  buckling  loads  were  obtained  by 
using  linear  anisotropic  shell  theory.  The  results  of  the  compres¬ 
sion  tests  indicate  that  the  cylinders  buckled  at  65  to  857o  of 
the  loads  predicted  by  classical  linear  analysis.  The  buckling 
mode  was  coupled  with  catastrophic  shear  failure  in  the  smaller 
cylinders.  The  larger  cylinders  buckled  into  the  classical 
diamond -shaped  pattern  with  the  basic  integrity  of  the  cylinder 
maintained , 
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The  composite  moduli  required  for  the  buckling  analysis  are 
determined  for  each  cylinder  from  experimental  test  data.  The 
data  are  obtained  by  the  classical  method  of  using  three  tests: 
torsion,  internal  pressure,  and  axial  tension  or  compression.  It 
is  demonstrated  analytically  that  the  anisotropic  coupling  be¬ 
tween  stretching  and  shear  only  influences  the  determination  of 
the  moduli  by  terms  of  the  order  of  the  square  of  the  ratio  of 
the  thickness  of  the  cylinder  to  the  radius  of  the  cylinder. 

Thus,  because  of  the  thinness  of  the  cylinders,  the  anisotropic 
coupling  has  a  negligible  effect  on  the  experimental  determina¬ 
tion  of  moduli.  Also,  the  moduli  obtained  from  the  experimental 

data  are  compared  with  those  predicted  by  the  theory  of  S.  Tsai.^ 
In  general,  agreement  is  good  when  a  value  of  0.3  is  assigned  to 
his  continguity  factor. 


II.  ANALYTICAL  PROCEDURE 


The  buckling  loads  for  clamped  filament -wound  cylinders  under 
uniform  axial  compression  are  analytically  determined  by  using 
linear  anisotropic  shell  theory.  The  general  procedure  for  ob¬ 
taining  buckling  solutions  for  anisotropic  cylinders  with  speci- 
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fied  boundary  conditions  has  been  given  by  Cheng  and  Ho.  For 
the  sake  of  completeness,  the  derivation  of  the  solution  is  given 
in  Appendix  B  and  is  briefly  outlined  here. 

The  displacement  equations  of  equilibrium  that  are  used  are 

obtained  by  combining  Flugge*s^^  stress  equations  of  equilibrium 
and  strain -displacement  expressions  with  an  anisotropic  stress - 
strain  relationship.  The  eighth  order  system  of  displacement 
equations  governing  the  stability  of  the  cylinder  are  solved  by 
first  assuming  a  wave  form  for  the  displacement  field  in  terms 

9  and  1 1 

of  the  axial  and  angular  coordinates.  The  requirement 

that  the  equations  of  equilibrium  be  satisfied  yields  a  character¬ 
istic  equation  governing  the  eight  axial  wave  numbers  that  are 
permitted  for  a  given  load.  Since  each  axial  wave  number  satisfies 
the  characteristic  equation,  the  displacement  modes  corresponding 
to  each  wave  number  can  be  superimposed.  The  modal  amplitudes  are 
governed  by  the  requirement  that  the  prescribed  end  conditions  be 
satisfied.  The  lowest  axial  load  that  simultaneously  satisfies 
the  differential  equations  of  equilibrium  and  the  end  conditions 
is  the  critical  buckling  load. 

The  numerical  solution  for  the  buckling  load  is  obtained  in 
the  following  manner.  First,  a  lower  bound  for  the  load  is  ob¬ 
tained  from  the  characteristic  equation  by  accepting  whatever 
boundary  conditions  occur  with  an  integral  number  of  half-waves 
specified  in  the  axial  direction.*  For  a  clamped  cylinder 
(v  =  V  =  u  -  Sw/5x  -  0  at  X  -  ±L/2),  the  buckling  load  is  raised. 
With  a  given  number  of  circumferential  waves  and  with  specified 
geometric  and  material  properties  of  the  clamped  cylinder,  the 


*For  the  particular  case  of  an  orthotropic  cylinder,  the  end 
conditions  for  the  lower  bound  solution  are  the  usual  mixed  or 
"simple'*  supports  on  a  right -circular  section,  that  is,  w  =  v  = 

N  =M  =0atx=  ±L/2.  With  a  general  anisotropic  medium, 

XX  XX 

the  end  conditions  corresponding  to  the  lower  bound  solution 
depend  on  the  degree  of  anisotropy. 
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buckling  load  is  found  by  an  interval  halving  procedure.  The 
lower  bound  solution  is  used  for  the  first  estimate.  The  initial 
increment  is  also  obtained  from  the  solutions  of  the  character¬ 
istic  equation  by  choosing  an  increment  that  is  less  in  magnitude 
than  the  difference  between  the  buckling  load  for  the  fundamental 
mode  and  the  next  higher  mode.  With  each  trial  buckling  load,  an 
eighth  order  polynomial  equation,  obtained  by  expanding  the 
characteristic  determinant,  is  solved  for  the  eight  axial  wave 
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numbers  by  using  Muller ^s  method.  After  the  wave  numbers  are 
computed,  the  boundary  determinant  obtained  by  requiring  clamped 
conditions  is  evaluated  by  Grout  reduction.  The  iteration  is 
continued  until  the  sign  of  the  boundary  determinant  changes.  The 
interval  is  then  halved,  and  the  process  is  repeated  until  suffi¬ 
cient  accuracy  is  achieved.  Buckling  loads  corresponding  to 
different  numbers  of  circumferential  waves  can  be  found,  if  nec¬ 
essary,  to  determine  the  minimum  buckling  load. 

The  stiffnesses  required  for  the  buckling  analysis  are  ob¬ 
tained  from  three  simple  experiments.  The  tests  used  are  internal 
pressure,  axial  tension  or  axial  compression,  and  torsion,  with 
principal  strains  measured  throughout.  The  load-strain  slopes 
from  the  three  tests  provide  a  direct  determination  of  the  in¬ 
plane  stiffness  matrix.  By  definition,  the  in-plane  stiffness 
matrix  is  the  integral  of  the  composite  moduli  over  the  thickness 
of  the  shell.  The  integration  yields  four  simultaneous,  symmetric, 
algebraic,  nonhomogeneous  equations  relating  the  principal  compos¬ 
ite  orthotropic  moduli  for  f ilament -wound  construction  to  the  in¬ 
plane  stiffnesses.  The  equations  are  simply  inverted  by  Grout 
reduction.  The  moduli  thus  determined  are  used  to  compute  the 
bending  stiffness  and  the  coupling  stiffness  matrices  used  in  the 
buckling  analysis.  The  derivation  of  the  equations  required  to 
compute  the  stiffnesses  and  moduli  from  the  test  data  is  given 
in  Appendix  G. 


III.  EXPER IMENTAL  PROCEDURE 


The  general  procedure  used  for  determining  stiffnesses  and 
the  buckling  strength  of  each  cylinder  was  as  follows.  End 
plates  were  shrink  fitted  to  the  cylinder  and  external  ring 
clamping  was  applied.  The  ten  to  twelve  strain  gages  on  the 
cylinder  were  used  to  first  record  strains  up  to  500  micro  inches /in. 
in  each  of  the  three  tests  used  to  determine  initial  stiffnesses, 
and  then  were  used  to  record  strains  during  the  buckling  test. 

The  locations  and  orientations  of  the  strain  gages  are  given  in 
Fig.  1.  The  internal  pressure  test  was  conducted  first.  This 
test  provided  a  good  check  on  strain  gage  performance  as  well  as 
serving  its  intended  purpose  for  determining  stiffness.  Torsion 
arms  were  then  fitted  to  the  end  plates,  and  a  pure  torque  was 
applied  to  the  cylinder  in  a  torsion  frame.  The  torque  arms  were 
then  removed  and  the  cylinder  was  placed  in  a  universal  testing 
machine  for  the  axial  tension  test.  Afterward,  because  of  the ^ 
possibility  of  slippage  during  the  tension  test,  the  end  clamping 
rings  were  loosened  to  permit  reseating  the  end  plates  against 
the  ends  of  the  cylinder.  The  end  plates  were  then  reclamped 
and  the  cylinder  tested  in  axial  compression  until  buckling 
occurred . 

After  the  buckling  test,  four  photomicrographs  were  taken  of 
the  cross  section  of  the  cylinder  at  widely  spaced  locations.  The 
photomicrographs  provided  layer  thicknesses  and  an  estimate  of 
winding  quality.  If  buckling  occurred  in  a  confined  region,  the 
photomicrographs  were  taken  only  in  that  region.  The  geometric 
parameters  of  each  cylinder  are  listed  in  Table  1. 

Before  each  test,  all  strain  gage  channels  were  zeroed  with 
a  switching  and  balancing  unit.  After  each  increment  of  load, 
the  gages  were  run  through  quickly  to  prevent  gage  heating  and 
subsequent  zero  drift.  This  technique  was  followed  because  the 
filament -wound  composite  has  a  low  coefficient  of  thermal  con¬ 
ductivity.  The  technique  was  found  to  be  generally  successful 
in  preventing  zero  drift.  Raw  data  were  plotted  in  the  form  of 
pressure,  torque,  or  force  versus  strain  for  each  strain  gage 
read.  The  slopes  were  then  picked  from  the  graphs. 
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IV.  RESULTS 


A.  EXPERIMENTAL  RESULTS 


A  summary  of  the  averaged  load-strain  slopes  used  to  determine 
stiffnesses  and  moduli  for  each  cylinder  is  given  in  Table  2.  An 
indication  of  the  individual  results  for  each  strain  gage, 
throughout  the  complete  sequence  of  testing,  is  given  for  Cylinder 
9  (D  =  24  in.)  in  Fig.  2  thru  5.  The  nature  of  the  buckling 
failures  in  the  6,  12,  18,  and  24-in.  dia*  cylinders  is  indicated 
in  Fig.  6  thru  12. 

The  behavior  of  the  specimens  with  internal  pressure  loading 
was  as  expected.  All  strain  gages  indicated  extensional  deforma¬ 
tion  with  the  largest  output  from  the  circumferential  gages.  A 
netting  analysis  would  predict  equal  strains  for  the  longitudinal 
and  circumferential  directions.  The  isotropic  influence  of  the 
resin,  however,  produced  a  composite  result  between  a  balanced 
design  (equal  strain)  and  an  isotropic  material  (1:2  ratio  for 
strain).  The  strain  gages  generally  indicated  a  linear  behavior 
with  negligible  hysteresis.  Hysteresis  was  quite  noticeable, 
however,  in  Cylinders  8  and  9  (Fig.  2).  This  may  have  been  due 
to  not  letting  the  pressure  gage  stabilize  sufficiently  during 
unloading  for  the  tanks  with  large  volume,  or  it  may  have  been 
due  to  strain  gage  drift. 

In  general,  the  data  from  the  tension  stiffness  tests  indicate 
a  remarkably  linear  behavior  with  little  hysteresis.  In  this 
respect,  the  tension  data  of  Gages  4  and  5  of  Cylinder  9  (Fig.  3) 
are  not  typical  of  the  gages  on  the  other  cylinders.  The  results 
of  the  torsion  tests  also  indicate  little  hysteresis  or  non¬ 
linearity.  The  results  of  the  torsion  test  on  Cylinder  9  are 
shown  in  Fig.  4,  with  the  principal  compression  and  tension 
strains  plotted  in  the  same  direction  to  facilitate  a  comparison 
of  the  slopes.  It  should  be  noted  that  some  of  the  cylinders 
were  loaded  to  full  test  pressure,  tension  and/or  torsion  before 
the  actual  stiffness  test  was  conducted.  Some  cylinders  were  not 
preloaded.  However,  preloading  appeared  to  have  no  effect  on  the 
quality  of  the  results,  in  either  linearity  or  hysteresis. 


*Diame ter . 


With  the  larger  diameter  cylinders,  there  was  some  difficulty 
in  applying  the  tension  loads  required  to  achieve  500  micro  inches/ 
in.  in  the  axial  strain  gages  because  of  slippage  of  the  end 
plates.  The  clamping  force  was  then  considerably  increased  by 
additional  tightening  of  the  clamping  screws  so  that  no  slippage 
occurred.  Unfortunately,  keeping  the  same  degree  of  clamping 
during  the  compression  test  appeared  to  affect  the  buckling  load 
and  the  character  of  the  buckles.  This  is  explained  in  more 
detail  later. 

The  axial  compression  loads  that  buckled  the  cylinders  are 
listed  in  Table  2.  For  all  the  cylinders  except  6,  7,  and  8, 
the  values  are  as  expected.  The  failures  were  sudden  and  without 
warning;  loads  dropped  off  precipitously.  All  four  6 -in.  diameter 
cylinders  experienced  a  catastrophic  shear  failure  coupled  with 
diamond -shaped  buckles.  The  plastic  behavior  of  the  composite 
was  most  noticeable  in  Cylinder  1,  which  failed  at  14,020  lb 
while  the  load  was  being  maintained  for  the  recording  of  dial  and 
strain  gage  readings. The  6 -in.  dia  cylinders  with  nominal  L/D 
=  5  had  buckles  only  near  an  end  clamping  plate  (Fig.  6).  Figure 
7  shows  a  view  of  the  inside  of  the  same  cylinder  after  the  end 
plates  were  removed.  The  6 -in.  dia  cylinders  with  nominal  L/D  =  1 
had  buckles  that  covered  one -third  to  one -half  the  length  of  the 
cylinders.  Figure  8  is  a  view  of  Cylinder  1,  unloaded,  with  the 
buckle  diamonds  outlined  with  a  china -marking  pencil. 

In  the  larger  cylinders,  the  diamond -shaped  buckles  occurred 
with  little  visual  evidence  of  permanent  damage.  Figure  9  is  a 
view  of  Cylinder  4,  just  after  buckling,  showing  general  absence 
of  wall  rupture. t  In  no  case  did  the  diamond  buckles  occur  com¬ 
pletely  around  the  cylinder.  Generally,  they  only  encompassed 
about  three-quarters  of  the  circumference.  This  is  due,  of 
course,  to  slight  variations  in  the  wall  construction  that  trigger 
the  buckle  and  to  the  use  of  a  swivel  loading  head.  Once  buckling 
starts,  the  loading  head  twists  and  tends  to  relieve  the  unbuckled 
portions.  Within  groups  of  cylinders  of  the  same  size,  buckling 


^Strain  gage  readings  were  taken  after  every  1000 -lb  increment 
of  compression  load  because  of  the  usefulness  of  the  gages  in  de¬ 
termining  compression  moduli.  After  the  experience  with  Cylinder 
1,  however,  this  practice  was  followed  only  up  to  9000  to  13,000 
lb,  and  then  the  test  was  run  continuously  to  failure  with  only 
one  gage  monitored. 

tThe  shear  effects  were  absent  because  the  stresses  at 
buckling  were  halved  by  doubling  the  diameter. 


always  occurred  at  the  same  end  (except  for  Cylinder  7,  which  was 
influenced  by  the  end  clamping).  At  first  this  was  thought  to 
be  a  function  of  the  orientation  of  the  specimens  in  the  testing 
machine,  since  the  failures  occurred  at  the  end  loaded  by  the 
swivel  head.  But  changes  in  orientation  resulted  in  continued 
failure  at  the  same  end  with  respect  to  the  original  location  on 
the  winding  mandrel.  Thus,  it  was  concluded  that  the  location 
of  failure  was  a  function  of  the  winding  technique  and  orienta¬ 
tion  . 

The  18 -in,  dia  cylinders  exhibited  a  substantially  lower 
buckling  strength  and  a  mode  shape  that  occurred  in  a  narrow 
region  around  one  clamping  plate  (Fig.  10).  Buckling  was  not 
sudden  and  occurred  almost  inaudibly.  The  second  18 -in.  dia 
cylinder  was  instrumented  with  two  strain  gages  near  a  clamped 
end;  the  gages  turned  out  to  be  in  a  region  that  did  not  expe¬ 
rience  the  local  end  buckling.  The  two  gages  not  only  indicated 
a  considerably  different  slope  than  the  other  axial  gages,  they 
also  recorded  a  permanent  strain.  The  loading  of  the  18-in.  dia 
cylinders  was  discontinued  after  the  first  local  buckles  occurred. 
The  damaged  ends  were  cut  off  and  the  rest  of  the  two  cylinders 
were  saved  for  further  testing  (the  testing  of  the  shortened 
cylinders  is  described  later). 

The  first  24-in.  dia  cylinder  (Cylinder  8)  experienced  the 
same  local  end  buckling  result.  The  load  was  then  continuously 
increased,  however,  to  find  out  what  would  happen.  The  initial 
end  buckles  grew  continuously  larger  until  an  audible  snap 
occurred  and  the  load  dropped.  Table  2  shows  the  two  values  of 
buckling  load  for  Cylinder  8  when  buckling  was  first  noticed, 
and  when  the  configuration  suddenly  changed  and  the  load  dropped 
off.  In  Fig,  11  the  growth  of  the  buckle  diamonds  is  shown  by 
the  letters  in  china -marking  pencil.  The  first  buckle  occurred 
at  "A."  The  state  corresponding  to  the  higher  load  in  Table  2 
is  labeled  ’’E."  Upon  unloading,  the  buckles  disappeared  gradu¬ 
ally,  rather  than  snapping  out. 

The  local  end  buckling  phenomenon  was  first  suspected  to  be 
caused  by  a  lower  stiffness  near  the  ends  (photomicrographs  taken 
afterward  in  the  buckled  region,  however,  gave  no  support  to  this 
suspicion).  Then,  an  inspection  of  the  photographs  of  the  end 
buckles  showed  they  occurred  directly  above  the  clamping  screws 
(Fig.  10).  The  consistency  of  this  one-to-one  correspondence  in¬ 
dicated  that  the  extreme  degree  of  clamping  warped  the  cylinder 
out  of  its  initial  shape  and  created  an  initially  imperfect  geome¬ 
try  that  seriously  influenced  the  specimen's  buckling  resistance. 
The  split  clamping  rings  for  the  large  diameter  cylinders  were 
apparently  too  flexible  to  distribute  the  clamping  pressure  uni¬ 
formly  . 


The  clamping  procedure  was  corrected  before  the  second  24-in. 
dia  cylinder  (No.  9)  was  tested.  The  number  of  clamping  screws 
was  doubled  to  24  and  the  screws  were  not  made  as  tight.  The 
cylinder  reached  its  proper  buckling  load.  No  local  end  buckles 
occurred.  The  failure  was  sudden  and  the  cylinder  buckled  into 
the  classical  diamond  pattern  (Fig.  12).  The  buckles  also  dis¬ 
appeared  suddenly  during  unloading.  Some  of  the  extreme  corners 
of  the  diamond  buckles  exhibited  layer  separation.  A  photomicro¬ 
graph  taken  of  the  cylinder  cross  section  is  shown  in  Fig.  13. 

In  addition  to  providing  layer  thickness  data  and  a  qualitative 
description  of  the  caliber  of  construction,  the  photomicrograph 
also  indicated  a  crack  that  was  not  visually  observed  after  the 
test.  Figure  14  provides  a  detailed  magnification  of  the  crack. 

The  last  series  of  tests  were  on  the  two  18 -in.  dia  cylinders 
whose  damaged  ends  were  cut  off.  The  new  nominal  L/D  ratio  was 
0.76.  The  cylinder  with  cp  =  0  was  somewhat  damaged  during  the 
cutting  process.  Despite  this,  the  cylinders  were  subjected  to 
the  complete  sequence  of  testing,  and  the  three  stiffness  tests 
indicated  little  change  in  stiffness.  The  clamping  procedure 
applied  to  the  shortened  18 -in.  dia  cylinders  was  similar  to  that 
on  Cylinder  9.  The  number  of  clamping  screws  was  doubled  to  24 
and  the  screws  were  not  made  as  tight.  The  buckling  strengths 
of  the  two  18 -in.  dia  cylinders  were  found  to  be  51  and  897o 
greater  than  the  strengths  of  the  same  cylinders  using  the  previous 
clamping  procedure.'^  No  end  buckles  occurred  during  the  axial 
compression  tests.  The  failures  were  sudden  and  the  cylinders 
buckled  into  the  classical  diamond  pattern  (Fig.  15). 


B.  ANALYTICAL  RESULTS 


The  analysis  of  the  data  obtained  in  this  study  is  summarized 
in  Tables  3  thru  6. 


*With  a  cylinder  whose  length  is  of  the  order  of  magnitude  of 
the  diameter,  small  changes  in  length  have  negligible  influence 
on  buckling  load. 
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The  moduli  obtained  from  the  experimental  data  are  given  in 

Table  3.  First,  the  in-plane  shear  stiffness  was  obtained 

66 

from  an  average  of  the  principal  tensile  and  principal  compressive 
small  strain  data  in  the  torsion  test.  Next,  the  in -plane  stiff¬ 
nesses  ^22^  ^11’  ^12  calculated  by  combining  the  small 

tensile  strain  data  from  the  internal  pressure  test  with  the  small 
strain  data  from  the  axial  tension  test,  the  small  strain  data 
from  the  axial  compression  test,  and  the  higher  strain  data  from 
the  axial  compression  test.  Finally,  moduli  were  calculated  from 
Eq  [C-13]  and  [C-14]  in  Appendix  C.  The  table  theoretically  in¬ 
dicates  a  dependence  of  the  shear  modulus  on  the  axial  tension  or 
compression  tests  only  because  of  the  slight  coupling  when  cp  ^  0 
(the  anisotropic  cylinders). 

The  equations  used  for  determining  moduli  from  stiffnesses 
for  cp  =  0  (the  orthotropic  cylinders)  offer  direct  insight  into 
which  moduli  can  be  computed  accurately.  For  Cp  =  0  the  equations 
are : 

<=22'  ■  -  VllVlC^’-c)'  -(’’p/l 

'=11°  ■  (''“c-'ll  -  V22)/((2»c)'  -  ("p)'] 

"n'  -  "ll/" 


<=66°  - 

where  h  is  the  thickness  of  a  circumferential  wrap  layer,  h  is 
c  p 

the  thickness  of  the  axial  wrap  layer,  and  h  is  the  total  thick¬ 
ness  of  the  cylinder  wall. 

c  c  c 

The  moduli  ^  ,  and  are  obviously  insensitive 

22  12  66 

to  slight  errors  in  the  stiffnesses.  However,  the  modulus 

Q 

perpendicular  to  the  fibers,  ,  is  small  and  is  determined  by 

the  difference  between  two  large  numbers.  This  is  the  reason  for 
the  greatest  fluctuation  in  data  in  Table  3,  appearing  in  the 
Young *s  modulus  perpendicular  to  the  fibers.  The  data  for  shear 
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modulus  show  the  greatest  consistency  of  all  the  moduli,  because 
they  are  obtained  from  a  direct  calculation.  Young ^s  modulus  in 
the  direction  of  the  fibers  is  not  quite  as  consistently  obtained, 
but  the  fluctuations  are  tolerable.  The  data  for  Poisson's  ratios 
can  only  be  understood  by  recalling  the  definitions: 


Major  Poisson's  ratio  =  ^ 

Minor  Poisson's  ratio  =  * 


The  great  fluctuations  indicated  in  the  major  Poisson's  ratio 
results  are  caused  by  the  fluctuations  in  the  modulus  perpendicular 

to  the  fibers.  The  data  for  remarkably  consistent  (the 

Q 

degree  of  consistency  of  indicated  by  the  consistency  of 

the  minor  Poisson's  ratio  column).  Since  is  the  quantity  of 

interest  in  a  shell  analysis,  the  effect  of  the  fluctuations  in 
the  major  Poisson's  ratio  data  is  more  illusory  than  real. 

A  survey  of  the  results  in  Table  3  indicates  no  consistent 
ratio  of  the  moduli  obtained  by  using  the  initial  axial  tension 
data  and  the  moduli  obtained  by  using  the  initial  axial  compres- 
s ion  data.  However,  the  differences  obtained  by  using  the  two 
sets  of  initial  axial  test  data  are  small.  Of  course,  greater 
differences  in  moduli  are  obtained  by  using  the  higher  strain 
data  of  the  axial  compression  test. 

Analytical  predictions  of  the  cylinder  buckling  loads  and  a 
comparison  with  the  experimentally  determined  loads  are  presented 
in  Table  4.  The  analytical  predictions  are  based  on  stiffnesses 
obtained  by  using  the  high  strain  axial  compression  data  combined 
with  the  data  from  the  internal  pressure  and  torsion  tests.  In 
general,  the  buckling  predictions  from  the  characteristic  equation 
(corresponding  to  a  simple  support  solution  when  cp  =  0)  are  negli¬ 
gibly  different  from  the  predictions  of  shell  analysis  with  clamped 
ends  specified.  The  magnitudes  of  the  wave  numbers  encountered 
make  this  physically  obvious. 

When  classical  snap  buckling  occurred  the  experimental  loads 
were  65  to  857o  of  the  linear  theory  predictions.  When  local  end 
buckles  occurred,  due  to  nonuniform  clamping  that  warped  the  ini¬ 
tial  end  shape  of  the  cylinders,  the  experimental  loads  were 


40  to  43%  of  the  classical  linear  theory  predictions.  The  tabu¬ 
lations  given  in  Table  4  for  Cylinders  6  and  7  correspond  to  the 
original  18 -in.  dia  cylinders  that  experienced  local  end  buckles. 
The  results  given  for  Cylinders  6a  and  7a  correspond  to  the  same 
18 -in.  dia  cylinders  with  their  damaged  ends  cut  off.  Cylinders 
6a  and  7a  were  tested  with  the  improved  clamping  procedure.  The 
two  ratios  tabulated  for  Cylinder  8  correspond  to  the  first 
occurrence  of  the  end  buckles  and  to  the  final  snap  after  more 
load  was  applied.  The  two  6 -in.  dia  cylinders  with  L/D  ~  5 
behaved  in  a  manner  similar  to  the  6 -in.  dia  cylinders  with  L/D 
~  1.  For  the  two  long  cylinders,  the  buckling  load  corresponding 
to  n  =  2,  m  =  1  was  just  slightly  higher  than  the  load  correspond¬ 
ing  to  n  =  8,  m  ~  40.  Cylinder  2  was  a  6 -in.  dia  cylinder  that 
was  initially  rejected  because  of  noticeable  wetting  in  the  com¬ 
posite  construction.  The  wetting  appeared  to  have  little  effect, 
however,  on  the  buckling  strength  of  the  cylinder.  Before  the 
cylinders  were  tested,  the  larger  diameter  cylinders  were  expected 
to  show  a  poorer  comparison  with  linear  theory  than  would  the 
smaller  cylinders.  This  was  based  on  the  increased  sensitivity  of 
the  larger  diameter  cylinders  to  initial  imperfections.  The 
larger  cylinders  sustained  remarkably  high  loads,  however,  and 
generally  exhibited  the  same  percentage  of  classical  load  as  the 
smaller  ones.  This  was  due  to  winding  the  cylinders  around  ac¬ 
curately  machined  aluminum  mandrels  and,  presumably,  also  was  due 
to  the  circumferential  stiffening. 

The  last  column  in  Table  4  indicates  the  ratio  of  the  buckling 
load  predicted  by  orthotropic  theory  (with  coupling  between 
stretching  and  shear  ignored)  to  the  buckling  load  predicted  by 
anisotropic  theory.  The  effect  of  coupling  between  stretching 
and  shear  is  obviously  negligible  for  the  cylinders  fabricated  in 
this  study.  Indeed,  the  cylinders  were  purposely  designed  for 
small  coupling.  The  and  B^^  elements  of  the  [B]  matrix, 

created  by  the  polar  winding,  actually  alternate  in  sign  at  each 
of  two  weave  lines  in  a  plane  that  is  skew  with  respect  to  the 
longitudinal  axis  by  the  angle  Cp  (Appendix  C).  If  the  coupling 
between  stretching  and  shear  were  large,  the  nonhomogeneity  of  the 
[B]  matrix  would  make  the  analysis  of  Cheng  and  Ho  theoretically 
inapplicable.  The  theoretical  difficulty  is  probably  not 
important,  however,  for  the  following  reasons: 

1)  If  B^^  and  B^^  were  homogeneous,  an  arbitrary  change 

in  their  sign  would  not  affect  the  buckling  load  of 
a  cylinder  with  homogeneous  boundary  conditions  speci¬ 
fied  on  right -circular  ends; 
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2)  Buckling  is  of  a  local  character  because  of  the  high 
number  of  circumferential  and  axial  waves  in  the  buck¬ 
ling  mode. 

For  these  two  reasons,  two  diametrically  opposite  regions  on  the 
cylinder  (that  are  not  at  the  weave  lines)  will  buckle  at  the  same 
applied  load.  If  an  analysis  were  to  be  made  that  would  include 
the  nonhomogeneity  of  the  [B]  matrix,  it  would  probably  show  a 
slight  effect  from  a  rapidly  decaying  boundary  layer  near  the  two 
weave  lines. 

The  effect  of  using  initial  tension  moduli  to  calculate  buckling 
loads,  instead  of  the  high  strain  axial  compression  data,  is  indi¬ 
cated  in  Table  5.  Generally,  the  difference  between  the  two  sets 
of  calculations  is  not  large.  The  largest  effect  occurs  in  Cylinders 
3  and  4.  The  10  to  14%  effect  on  Cylinders  3  and  4  is  predominantly 
due  to  the  difference  between  the  moduli  perpendicular  to  the  fibers. 
Although  the  modulus  perpendicular  to  the  fibers  is  small,  it  has 
an  important  influence  on  D-.  ^  ,  which  in  turn  has  an  important  in- 

4 

fluence  on  the  buckling  strength  of  the  cylinders  . 

There  are  many  analytical  models  available  for  the  prediction 
of  elastic  moduli  of  unidirectional  fiber -re inforced  composites. 

The  selection  of  the  proper  model  is  subject  to  the  constraint  im¬ 
posed  by  the  geometric  arrangem.ent  of  the  glass  fibers  within  the 
resin  matrix.  Each  of  the  layers  in  the  cylinders  fabricated  for 
this  buckling  study  generally  have  fibers  that  are  randomly  distri¬ 
buted.  In  addition,  many  of  the  fibers  are  contiguous  because  of 
the  low  resin  content.  For  these  reasons,  the  analytical  model  of 
8 

S.  Tsai  appears  to  be  an  appropriate  model  to  compare  with  the 
moduli  based  on  the  experimental  data.  Young's  modulus  in  the  di¬ 
rection  of  fibers  is  based  on  the  usual  assumption  that  the  fibers 
and  the  resin  carry  loads  proportional  to  their  stiffnesses.  The 
models  for  Young’s  modulus  perpendicular  to  the  fibers,  major 
Poisson's  ratio,  and  in-plane  shear  modulus  are  constructed  by  linear 
interpolation  between  two  extreme  cases:  (1)  all  fibers  are  isolated 

by  resin;  (2)  all  fibers  are  contiguous.  S.  Tsai's  contiguity  fac¬ 
tor  (C)  indicates  the  degree  to  which  the  fibers  touch.  The  case 
C  =  0  indicates  all  fibers  are  isolated.  The  case  C  =  1  indicates 
all  fibers  are  contiguous. 


Table  6  provides  a  comparison  of  the  experimentally  determined 
initial  tension  moduli  with  the  analytical  model.  The  parameters 
used  in  the  analytical  computations  are 

=  12.4  X  10^  psi 

=  0.22 

7f  =  2.49 

E  =  0 . 5  X  10^  psi 
m 

V  =  0.35 
m 

7  =  1.22 

m 

k  =  1.0, 

where  E  is  Young's  modulus;  v  is  Poisson's  ratio;  7  is  specific 
gravity;  subscripts  f  and  m  refer  to  the  filament  and  matrix,  re¬ 
spectively;  and  k  is  a  filament  misalignment  factor.  The  con¬ 
tiguity  of  the  fibers  is  quite  noticeable  in  the  photomicrographs 
that  were  taken  and  has  a  noticeable  influence  on  Young's  modulus 
perpendicular  to  the  fibers  and  on  the  shear  modulus.  C  =  0.3  ap¬ 
pears  to  be  representative  of  the  cylinders  tested  in  this  study. 

The  contiguity  factor,  crude  as  it  is,  does  provide  a  simple  first 
approximation  for  the  effect  of  contiguity  of  the  fibers.* 

The  percent  resin  content  in  the  composite,  not  the  total  per¬ 
cent  resin  content,  is  used  in  the  theoretical  model  computations. 

The  total  resin  contents  were  determined  by  loss  of  weight  due  to 
ignition  of  samples  of  the  specimens  taken  after  the  tests  were 
completed.  The  values  of  resin  content  in  the  composite  material 
(Table  6)  were  determined  by  using  the  photomicrographs  to  estimate 
how  much  the  outer  resin  layer  contributed  to  the  total  resin  con¬ 
tent.  While  the  difference  between  the  two  sets  of  resin  contents 
is  small,  its  influence  is  felt  in  the  theoretical  shear  modulus. 

If  the  total  percent  resin  content  had  been  used,  the  contiguity 
factor  C  =  0.4  would  have  yielded  a  good  comparison  with  experi¬ 
mental  data,  rather  than  C  =  0.3.  This  sensitivity  of  shear  modulus 
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to  contiguity  has  been  noted  by  Ekvall  and  can  be  seen  in  S.  Tsai's 
results.  It  is  an  example  of  where  a  simple  test  can  provide  ac¬ 
curate  values  of  modulus  while  an  accurate  theoretical  prediction 
is  difficult  to  achieve, 

*For  filament  wound  composites  whose  fibers  are  isolated  by  the 

13 

resin  matrix,  the  method  of  Hashin  and  Rosen  offers  another  ra¬ 
tional  method  of  analysis.  In  addition  to  using  the  same  expression 
for  the  modulus  in  the  direction  of  fibers,  their  shear  modulus  is 
identically  equal  to  that  of  S.  Tsai  with  C  =  0. 
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V.  CONCLUDING  REMARKS 


For  cylinders  with  diameter  to  thickness  ratios  of  167  to 
643,  classical  linear  shell  buckling  theory  has  been  found  to 
provide  a  reasonable  basis  of  comparison  with  the  axial  compres¬ 
sion  loads  that  a  filament -wound  cylinder  can  sustain.  Still, 
the  solutions  to  the  following  problems  should  be  obtained: 

1)  A  large  deflection  analysis  for  anisotropic  cylinders 
that  includes  initial  geometric  imperfections.  When 
design  criteria  are  used  in  which  the  avoidance  of 
buckling  is  at  least  a  part,  the  layer  arrangements 
in  a  filament -wound  cylinder  should  be  chosen  to 
reduce  the  effect  of  initial  imperfections.  The 
large  deflection  analysis  that  includes  the  effect 

of  initial  imperfections  can  then  be  used  as  a  design 
tool ; 

2)  A  large  deflection  analysis  for  anisotropic  cylinders 
that  includes  the  effect  of  boundary  conditions  on 

the  prebuckling  stress  state.  Almroth^^  has  shown, 
for  isotropic  cylinders,  that  the  classical  buckling 
load  is  somewhat  reduced  when  the  effect  of  uniform 
clamping  on  the  prebuckling  state  is  considered; 

3)  Refined  methods  of  analysis  for  predicting  composite 
moduli.  The  effect  on  composite  moduli  of  fabrica¬ 
tion  processes,  a  detailed  study  of  filament  contig¬ 
uity,  and  large  strain  response  of  the  resin  matrix 
should  be  obtained. 
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*Total  length  between  clamping  rings  is 


Internal  Pressure  Torsion 
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Plus 


Axial  Test 
Slopes  Used 

Young's  Modulus 
in  Direction 
of  Fibers,  10^  psi 

Initial  Tension 

8.86 

Initial  Compres ■ 
s  ion 

8.64 

Final  Compres¬ 
sion 

8.85 

Initial  Tension 

10.13 

Initial  Compres- 
s  ion 

10.37 

Final  Compres¬ 
sion 

10.30 

initial  Tension 

10.00 

Initial  Compres** 
'  s  ion 
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Final  Compres¬ 
sion 
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Initial  Tension 

9.10 

Initial  Compres¬ 
sion 

8.85 

Final  Compres¬ 
sion 
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Initial  Tension 
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Initial  Compres¬ 
sion 

8.52 

Final  Compres¬ 
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Initial  Tension 
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Final  Compres¬ 
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Initial  Tension 

8,61 
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sion 

8.69 

Final  Compres¬ 

sion 

9.11 

of  axial  test  slope  for  orthotrop 

Table  4  Cylinder  Buckling  Loads 
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Table  5  Effect  of  Method  Used  for  Determining  Moduli 
on  Theoretical  Buckling  Loads 


Table  6  Comparison  of  Tension  Moduli  With  an  Analytical  Model 


Method 

Used 

Gylinder 

No. 

Total  Resin 
Gontent 
(Percent  by 
Weight) 

Resin  Content 
in  Composite 
(Percent 
by  Weight) 

Young's  Modulus 
in  Direction  of 

Fibers^  10  psi 

hi" 

Young's  Modulus 
Perpendicular  to 

Fibers^  10  psi 

h2 

Major 

Poisson's 

Ratio 

hb" 

Shear 

Modulus, 

10^  psi 

Experiment 

1 

20.4 

18.2 

8.85 

3.09 

0.32 

1.39 

Tsai  C  =  0.4 

C  =  0.3 

C  =  0.2 

C  =  0.0 

8.68 

3.05 

2.75 

2.45 

1.86 

0.24 

0.25 

0.25 

0.26 

1.62 

1.43 

1.23 

0.84 

Experiment 

2 

20.4 

17.0 

8.61 

4.24 

0.22 

1.49 

Tsai  C  =  0.4 

C  =  0.3 

C  =  0.2 

C  =  0.0 

8.90 

3.21 

2.90 

2.59 

1.97 

0.24 

0.25 

0.25 

0.25 

1.69 

1.49 

1.29 

0.89 

Experiment 

3 

18.7 

16.2 

9.79 

3.16 

0.32 

1.52 

Tsai  C  =  0.4 

C  =  0.3 

C  =  0.2 

C  =  0.0 

9.04 

3.32 

3.00 

2.68 

2.04 

0.24 

0.24 

0.25 

0.25 

1.74 

1.54 

1.34 

0.93 

Experiment 

4 

19.3 

18.7 

9.10 

3.59 

0.28 

1.47 

Tsai  C  =  0.4 

C  =  0.3 

C  -  0.2 

C  =  0.0 

8.60 

2.99 

2.70 

2.41 

1.82 

0.25 

0.25 

0.25 

0.26 

1.60 

1.40 

1.21 

0.82 

Experiment 

5 

19.3 

17.4 

10.13 

2.77 

0.41 

1.46 

Tsai  C  =  0.4 

C  =  0.3 

C  =  0.2 

C  =  0.0 

8.83 

3.16 

2.85 

2.54 

1.93 

0.24 

0.25 

0.25 

0.25 

1.67 

1.47 

1.27 

0.88 

Experiment 

6 

6a 

7 

7a 

21.3 

19.0 

9.34 

8.80 

8.33 

8.88 

3.02 

3.14 

3.22 

3.20 

0.34 

0.32 

0.26 

0.32 

1.41 

1.37 

1.32 

1.34 

Tsai  C  =  0.4 

C  =  0.3 

C  =  0.2 

C  =  0.0 

8.55 

2.95 

2.66 

2.37 

1.80 

0.25 

0.25 

0.25 

0.26 

1.58 

1.39 

1.20 

0.81 

Experiment 

8 

21.3 

t 

8.13 

3.34 

0.26 

1.54 

9 

20.4 

t 

9.42 

3.06 

0.29 

1.29 

Experiment 

11 

19.4 

16.6 

8.38 

3.76 

0.22 

1.56 

Tsai  C  =  0.4 

C  =  0.3 

C  =  0.2 

C  =  0.0 

8.96 

3.25 

2.94 

2.63 

2.00 

0.24 

0.25 

0.25 

0.25 

1.71 

1.51 

1.31 

0.91 

Experiment 

12 

19.4 

16.4 

8.50 

4.33 

0.29 

1.54 

Tsai  C  =  0.4 
C  =  0.3 

G  =  0.2 

G  =  0.0 

8.99 

3.28 

2.97 

2.65 

2.02 

0.24 

0.25 

0.25 

0.25 

1.73 

1.52 

1.32 

0.92 
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Cylinder  L  (in.)  A  (deg) 


Fig.  1  Strain  Gage  Locations 


Fig,  2  Cylinder  9,  Internal  Pressure  vs  Strain 


100 


Strain  (micro-in. / in. ) 


Fig.  3  Cylinder  9,  Tension  Load  vs  Strain 


Fig.  5  Cylinder  9,  Compression  Load  vs  Strain 


wiiiiiiiiiiiiiiii 


Fig,  9  Cylinder  4  (12-in, dia)  in  Postbuckled  State 
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.dia)  in  Postbuckled  State 
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Fig.  11  Cylinder  8  (24-in.  dia)  after  Test 


Fig.  13  Cylinder  9,  Photomicrograph  of  Cylinder  Wall 
(Taken  at  an  angle  of  19°30'  with  respect 
to  the  axis  of  the  cylinder,  SOX) 


Fig.  14  Cylinder  9,  Detail  of  Crack  Shown  in  Fig.  13 


APPENDIX  A 


FABRICATION  OF  CYLINDERS 


A  single  lot  of  pre-impregnated  fiberglas  roving  was  obtained 
from  U.S.  Polymeric  Chemical  Co.  for  the  fabrication  of  all  the 
cylinders.  This  material  was  20  end,  S994  roving  with  an  HTS  finish, 
and  was  impregnated  with  an  epoxy  resin  system  designated  EF  787 
by  U.S.  Polymeric  Chemical  Co. 

All  cylinders  were  wrapped  with  a  roving  tension  of  10  lb,  and 
consisted  of  three  layers  of  material.  The  first  and  third  layers 
were  wrapped  circumferentially  and  contained  400  ends  per  inch  roving. 
The  middle  layer  parameters  were  systematically  varied  and  are  listed 
in  Table  1.  The  nominal  wall  thickness  of  all  cylindrical  specimens 
was  0.036  in.  All  cylinders  were  cured  in  the  same  oven,  using  a 
cam  that  resulted  in  the  following  cure  cycle:  6  hr  from  200  to 
350°F  plus  16  hr  at  350°F.  Following  cure,  the  specimens  were  cut 
to  the  length  shown  in  Table  1.  The  layer  thicknesses  in  Table  1 
are  averages  of  values  taken  from  four  photomicrographs  of  cross- 
sections  cut  from  the  cylinder  after  testing.  The  values  are 
exclusive  of  a  thin  resin  layer  that  exists  on  the  outer  surface 
of  each  of  the  cylinders  due  to  compaction  of  the  rovings. 
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APPENDIX  B 


SOLUTION  OF  EQUATIONS  OF  EQUILIBRIUM 


Appendix  B  describes  the  analysis  used  to  obtain  the  buckling 
strengths  of  clamped  filament -wound  cylinders  under  axial  compres¬ 
sion. 

Stress  Equations  of  Equilibrium  -  A  cylindrical  shell  of  con¬ 
stant  thickness  h  and  mean  radius  R  is  referred  to  a  modified 
cylindrical  system  of  coordinates  x,  6,  and  z  as  indicated  in  Fig. 
B-1;  X  is  measured  along  the  axis  of  the  shell,  0  along  the  circum¬ 
ference,  and  z  in  the  direction  perpendicular  to  the  middle  surface. 


The  linear  stress  equations  of  equilibrium  ^  in  terms  of  the 
notation  of  this  report  are: 


Sn 


1 

^  +  i  _ 

3x  R  Se  XX 


2  _2  ^,2 


Sx' 


R  M 


-v  2,2 

66  ^  ,  9  x0  o  u  _ 

R  Sx  R  Sxd^  ’ 


B-1] 


i  ]_  i  ^^0  „i  ^  , 

R  S©  Sx  „2  o0  R  c^x  XX  2  2 
R  ox 


N 


,  ‘'00  S  /  ,  \  , 


R 


N 


X  0  S 
R 


[B-2] 


r)^M 

XX  ,  I  ^  ^00  ,  1  ^  ^ex  1  ^  X0 

,2  2  ^.2  R  S0dx  R  S0Sx 

OX  R  o0 


N 


+  n"  ^  + 

R  XX  -X  z 

ox 


^00  /  ^U  1  ^  .  1.  5^W  \  ^  ^X0  S 

R  I  Sx  R  d0  R  ^^2  )  R 


Fig.  B-1  Cylindrical  Coordinate  System  for  Composite  Shell 
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where 


[B-4] 


The  superscript  i  refers  to  the  initial  stress  state;  stress  re¬ 
sultants  and  moments  without  superscripts  represent  the  additional 
stresses  that  occur  with  buckling;  u,  v,  and  w  are  the  displacements 
in  the  axial,  circumferential,  and  radial  directions,  respectively, 
which  arise  due  to  buckling.  The  cylinders  to  be  considered  are 
thin  (h/2R  «  1).  This  has  been  used  extensively  to  reduce  the 

number  of  initial  stress  terms  possible.  The  initial  stress  re¬ 
sultants  used  for  axial  loading  are: 


N' 


XX 


LOAD 

TtD 


N' 


QQ 


= 


Constitutive  Equations  -  The  anisotropic  stress -strain  relation¬ 
ship  for  a  cylinder  in  a  state  of  plane  stress  is 


0 

XX 

“hi 

h2 

he” 

'  e 

XX 

h2 

C22 

he 

• 

^xe  _ 

.he 

he 

'^6e_ 
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with  the  Flligge  strain  expressions  given  by 


o  au  ^ 

^x6  R  bo  Sx’  x0 


^  t  j:_  ^  ^v\ 

R  \  d0Sx  2R  ^0  2  5x/' 


The  constitutive  equations  for  the  cylinder  are  obtained  by 
substituting  Eq  [B-5]  and  [B-6]  into  [B-41 .  The  expressions  are 

9 

given  by  Cheng  and  Ho  and  are  repeated  here. 
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with  [A]  ,  [B]  ,  and  [D]  defined  by 


[A] 


[B] 


[D] 


dz 


zdz 


2, 

z  dz 


>  [B-9] 


For  the  cylinders  considered  in  this  study,  the  stiffness  matrices 
reduce  to  the  forms 


Ai 

^12 

0 

0 

0 

"Ai 

A2 

0 

[A]  = 

A^2 

^22 

0 

,  [B]  = 

0 

0 

®26 

,  [D]  = 

A2 

°22 

0 

0 

0 

A6_ 

_^6 

^26 

0 

0 

0 

See  Eq  [C-3],  [C-4],  and  [C-5]  in  Appendix  C. 


Characteristic  Equation  -  The  displacement  equations  of  equilib¬ 
rium  are  obtained  by  substituting  Eq  [B-8]  and  [B-7]  into  Eq  [B-1] 
thru  [B-3] .  The  substitution  of  displacements  of  the  form 


u  = 


V  = 


A  sin(  ^  +  n6 
R 


1„(: 


B  sin(  +  nO 
K 


w  =  C  cos 


/Ax 

(r 


) 

) 

) 


into  the  equations  of  equilibrium  yields  the  characteristic  deter¬ 
minant 


11 

A2 

"13 

12 

A2 

"23 

=  0 

[B-10] 

13 

^23 

"33 
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with 


and  n=05  1,  2,  3,  ,  .  .  . 
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If  A  =  mjt  R/L,  m  =  I5  2,  3  .  .  .  the  solutions  correspond  to 
boundary  conditions  that  are  determined  by  the  trigonometric  form 
of  the  assumed  solutions.  The  solutions  to  be  found  by  this  proce¬ 
dure  are  obtained  by  expanding  the  characteristic  determinant  into 
the  form 


(4) 


-b  P.N  +  P  -  0. 
2  XX  3 


[B-11] 


Equation  [B-11]  is  solved  by  the  New ton-Raphson  procedure 
The  first  approximation  for  is  given  by 

XX 


17 


N' 


XX 


hl^2- 


The  iteration  is  completed  when 


(n^  )  -  (nM 


(n^  ) 


<  1  X  10 


The  calculations  are  performed  by  an  IBM  1620  computer.  For  a  given 
value  of  n  the  axial  buckling  load  is  printed  for  m  =  1  through  a 
maximum  specified  value  of  m.  The  minimum  buckling  load  for  the 
given  value  of  n  is  also  printed.  The  procedure  is  then  repeated, 
with  n  increased  by  unity,  until  the  maximum  specified  value  for  n 
is  achieved.  Generally,  a  maximum  of  two  iterations  is  required 

for  each  .  The  case  n  =  0  sometimes  requires  five  iterations. 
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Effect  of  Boundary  Conditions  -  The  solution  for  arbitrary  bound¬ 
ary  conditions  uses  displacements  of  the  form: 


u 
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s,n[—  + 


•  f  4- 


cosi 


A,  X 

k 
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[B-12] 
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where  A,  represents  the  roots  of  the  eighth  order  polynomial  equa- 

K 

tion  obtained  by  expanding  the  characteristic  determinant  Eq  [B-10] 
in  powers  of  A,  Since  each  root  satisfies  the  characteristic  equa¬ 
tion,  the  displacements  corresponding  to  each  root  can  be  super¬ 
imposed  . 

Because  the  coefficients  of  the  polynomial  equation  are  real, 
a  complex  root  must  be  accompanied  by  a  root  that  is  its  complex 
conjugate.  Consider  the  two  roots 


\  =  hk  +  ^^ik 


\+l  hk  "  ’■^ik* 


Since  the  displacements  are  real,  we  have 


"k  = 


C  ,  +  ic., 

rk  ik 


C,  ,  =  C  ,  -  iC  . 

k+1  rk  ik 


The  expansion  of  the  k  and  k+1  components  of  w  in  Eq  [B-12]  in 
terms  of  its  real  and  imaginary  parts  yields  the  real  sum 


where 


”k  "k+1  -  (“rkhk  +  ^Akhk) 

+  (-2C^khk  +  ^AkAk)  ““ 

hk  "  (\kl)“"'’(Akl) 

Ak  ■  “=  (''rk  I)  (hk  I) 

Ak  ■  ('"tk  I)  (^ik  I) 

Ak  ■  (hk  !)“"*■  (hk  I) 
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The  u  and  v  displacements  are  treated  in  a  similar  manner.  The 
modal  amplitudes  and  are  found  in  terms  of  from  the  equa¬ 
tions  of  equilibrium.  Let 


then 


where 


\  ■  (4  +  ‘4)  ‘=k 
\  ■  (4  “ik)  "k 

u,  +  u,  ,  =  /zc  ,  -  2C.,  U^.  )  cos  n0 

k  k+1  \  rk  cr  ik  ci/ 

+  he  ,  -  2C.,U^.)  sin  ng 

\  rk  sr  ik  si/ 
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cr 


SI 


f  t  -ft 

rk^3k  ik^4k 

f"^  t  +  f^  f 
ik^3k  rk^4k 

fA  .  A  u 

rk^Zk  ik^lk 

A  A 

f  t  -ft 

ik^2k  rk^lk* 


The  sum  +  ^k+1  obtained  by  replacing  superscript  A  by 

superscript  B.  Slopes  and  stress  resultants  are  obtained  by  suitable 

differentiations  of  the  displacements.  When  two  real  roots  occur, 

instead  of  a  pair  of  complex  conjugate  roots,  they  are  represented 

in  the  summations  for  the  total  deflections  by  the  coefficients  of 

2C  ,  with  2C..  =  0, 
rk  ik 


44 


All  of  the  ingredients  required  for  satisfying  boundary  condi¬ 
tions  are  now  available.  The  clamped  conditions,  corresponding  to 

the  experiments,  are  u  =  v  =  w  =  =  0  at  x  =  ±L/ 2 .  Setting  the 

coefficients  of  cos  ne  and  sin  n0  equal  to  zero,  one  arrives  at  eight 
simultaneous  homogeneous  algebraic  equations.  For  given  geometric 

and  stiffness  properties  of  the  cylinder,  the  lowest  value  of 

XX 

that  sets  the  determinant  of  the  coefficients  of  20^^  and  2C^^  equal 

to  zero  is  the  load  per  unit  of  circumferential  length  that  is  re¬ 
quired  to  buckle  the  cylinder. 
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APPENDIX  C 


DETERMINATION  OF  STIFFNESSES  AND  MODULI 


Appendix  C  describes  the  analysis  used  to  obtain  the  stiffnesses 
and  moduli  of  filament -wound  cylinders  from  suitable  test  data.  The 
tests  required  are  internal  pressure,  axial  tension  or  axial  com¬ 
pression,  and  torsion. 


Transformation  Equations  for  Composite  Moduli  -  The  composite 
modulus  matrix  for  a  filament -wound  structure  in  a  state  of  plane 
stress  has  an  orthotropic  form  when  referred  to  axes  along,  and 
perpendicular  to,  the  fiber  direction.  Let  x‘,  y',  z'  denote  such 
a  system  of  axes,  with  the  fibers  directed  along  x'.  The  ortho¬ 
tropic  modulus  matrix  C’  can  be  written  as 


[C’] 


pi  pi 

11  12 


pi  p  f 

12  22 


0  0 


0 
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An  in-plane  rotation  to  a  new  coordinate  system  (x,y,z)  is  ac¬ 
complished  by 
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cos  cp  ”sin  cp  0 

x' 
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i — \ 
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where  cp  is  the  angle  of  rotation  shown 


in  Fig.  C-1 . 
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Fig.  C-1  Transformation  Axes  for  Composite  Moduli 


Since  [C’]  transforms  according  to  the  standard  rules  for  a 
fourth  rank  tensor,  its  representation  in  the  new  coordinate  system 
is  given  by* 


where 


and 
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c  =  cos  cp 
s  =  sin  cp . 


Constitutive  Equations  -  The  stress -strain  relationship  for  a 
polar  half-layer  referred  to  the  longitudinal  and  circumferential 
directions  is 
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XX 

Cii  C^2 
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XX 
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[C-2] 


'"The  compact  representation  of  the  fourth  rank  tensor  components 

C!.,  T  by  matrix  notation  C’  is  only  a  convenience.  To  transform 
ijkl  mn 

the  components  of  the  moduli  to  other  axes  it  is  necessary  to  use 
the  complete  tensor  notation. 
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3^ 


where , 


t  h 

j  refers  to  the  j  polar  half-layer. 


sig.(z)  =  ±1  for  polar  angle  =  ±cp . 

The  wall  of  each  cylinder  fabricated  for  this  study  is  composed 
of  three  layers.  The  inner  and  outer  layers  are  circumferential 
windings  of  equal  thickness  h^5  the  middle  layer  is  a  polar  winding 

ith  two  half-layers  wound  at  angle  of  wrap  ±cp  and  total  thickness 
The  in-plane,  coupling,  and  bending  stiffness  matrices  [A], 

P’ 

[B],  and  [D]  are  obtained  by  applying  the  [C]  matrix  in  Eq  [C-2]  to 
Eq  [B-9] .  The  result  is 
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with 


It  is  assumed  that  the  quality  of  the  construction  of  the  cyl¬ 
inder  is  such  that  each  layer  has  the  same  uniform  properties  when 
referred  to  the  direction  of  winding.  This  can  be  qualitatively  con¬ 
firmed  for  an  individual  cylinder  by  using  a  light  source  to  illu¬ 
minate  the  wall  of  the  cylinder.  The  transparency  of  a  thin  fila¬ 
ment-wound  cylinder  permits  a  visual  inspection  of  the  uniformity 
of  construction.  The  lack  of  a  complete  bond  between  resin  and 
fibers  is  simply  detected  by  black  streaks  where  wetting  has  oc¬ 
curred  . 
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The  choice  of  the  sign  in  [B]  depends  on  which  side  of  the  polar 
wrap  is  considered.  While  the  sign  convention  chosen  for  cp  is  arbi¬ 
trary  ^hence  the  sign  convention  of  is  arbitrary^,  the  ex¬ 

istence  of  the  two  different  signs  is  not  arbitrary.  The  sign  in 
[B]  alternates  at  each  of  two  diametrically  opposite  weave  lines  in 
the  polar  wrap.  The  two  weave  lines  are  parallel  to  each  other  in 
a  plane  that  intersects  the  center  of  the  cylinder  and  is  skew  with 
respect  to  the  longitudinal  axis  by  the  angle  cp. 

In  order  to  understand  why  the  weave  occurs,  one  must  understand 
how  a  polar  wrap  is  wound  around  a  mandrel.  The  plane  in  which  the 
fibers  are  wound  is  kept  constant  while  the  mandrel  is  continuously 
rotated  until  a  complete  rotation  of  360  deg  is  achieved.  Figures 
C-2  thru  C-7  illustrate  three  phases  of  the  winding  sequence." 
Figures  C-2  and  C-3  are  front  and  rear  views,  respectively,  of  the 
first  polar  encirclement  of  roving  around  the  mandrel.  Notice  that 
when  the  roving  is  referenced  with  respect  to  the  local  surface  of 
the  mandrel  two  diametrically  opposite  fibers  in  the  roving  are 
oriented  at  angles  that  are  equal  to  each  other  in  magnitude  but  are 
opposite  in  sign.  Figures  C-4  and  C-5  represent  front  and  reart 
views  of  the  rovings  on  the  mandrel  after  the  mandrel  has  rotated 
180  deg.  One  can  see  the  beginning  of  the  overlapping  of  rovings. 
Also,  by  considering  only  those  parts  of  the  rovings  that  are  in 
contact  with  the  surface  of  the  mandrel,  one  can  see  that  the  half- 
layer  on  the  left  side  of  Fig.  C-4  has  an  orientation  that  is  op¬ 
posite  to  the  orientation  of  the  half-layer  on  the  right  side  of 
Fig.  C-4.  The  same  is  true  for  the  rear  view  shown  in  Fig.  C-5. 
Figures  C-6  and  C-7  are  front  and  rear  views  of  roving  on  the  man¬ 
drel  just  before  the  mandrel  has  completed  its  360  deg  rotation 
(hence,  just  before  the  two  outside  half-layers  and  the  two  weave 
lines  are  completed).  A  detailed  view  of  a  weave  "line*'  in 
Cylinder  1  is  given  by  the  photomicrograph  in  Fig.  C-8.  The 
rovings  in  the  outer  polar  half  layer  are  seen  to  gradually  turn 
into  what  is  called  the  inner  polar  half  layer.  Also,  the 
transition  length  required  for  this  is  seen  to  be  of  the  order 
of  magnitude  of  the  thickness  of  the  shell. 


*A  black-coated  mandrel  and  gaps  between  rovings  are  used  in 
Fig.  C-2  thru  C-7  for  illustration  purposes.  The  cylinders  fabri¬ 
cated  for  the  buckling  study  were  wound  around  cylindrical  aluminum 
mandrels  and  the  rovings  were  laid  adjacent  to  each  other. 

tThe  front  views  in  Fig.  C-2,  C-4,  and  C-6  are  taken  from  the 
same  position.  Similarly,  the  rear  views  in  Fig.  C-3,  C-5,  and 
C-7  are  taken  from  the  same  rear  position,  directly  opposite  the 
front  view. 
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Determination  of  Stiffness  and  Moduli  -  The  simplified  forms 
of  the  stiffness  matrices  given  by  Eq  [0-3],  [C-4],  and  [C-5]  reduce 
the  constitutive  relations  Eq  [B-8]  to 
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ith  similar  expressions  for  and  the  bending  moments. 


For  the  moment,  assume  [B]  is  homogeneous.  Since  N  =  constant, 
N  =  constant,  and  N  ^  =  constant  in  the  three  tests  considered, 

66  xe 

there  can  be  no  changes  in  displacements  with  respect  to  0.  If  all 
strain  data  are  obtained  at  the  midsection  of  the  cylinder  the 
curvature  due  to  restraint  of  the  cylinder  ends  can  be  neglected 
and  thus: 
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Eq  [C-8]  reduces  to 
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Now  consider  longitudinal  tension  and  internal  pressure  applied 

to  the  shell.  With  no  restraint  provided  to  rotation  at  the  ends, 

N  ^  =  0,  Therefore, 
xG 
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The  expressions  in  brackets  differ  from  unity  by  terms  of  order 
2 

(h/R)  .  Because  of  the  thinness  of  the  shells,  the  difference  can 
be  effectively  neglected.  The  constitutive  equations  used  for 
determining  ^22  orthotropic  cylinders  (cp  =  0) 

are  therefore  equally  valid  for  the  anisotropic  cylinders  (cp  0)  , 
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The  in-plane  stiffnesses  k^^'^are  determined  by 

first  applying  uniform  tensile  loads  to  the  ends  of  the  cylinder 


N 

XX 

and 


=  Load/nD,  =  0  and  measuring  the  principal  strains  ) 

/  o  \  ^ 

y  em  *  '^ft;er  the  tension  test  is  completed,  internal  pressure 

p  is  applied  =  pR,  and  the  principal  strains 

and  1^00  I  3.re  measured.  The  tests  provide  four  equations  for 

determining  the  three  stiffnesses.  Three  of  the  equations  can  be 
combined  to  yield 


[C-10] 


The  fourth  equation  from  the  tests  provides  a  check  on 


[C^ll] 
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The  same  procedure  is  applied  to  the  torsion  test.  With  no 
axial  loads  or  internal  pressure,  Eq  [C-9]  yields 


Again,  for  (h/R)^  «  1,  the  bracket  term  is  effectively  unity  and 
the  constitutive  equation  used  for  determining  for  the  ortho¬ 

tropic  cylinders  is  equally  valid  for  the  anisotropic  cylinders, 
that  is, 
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The  engineering  shear  strain 


Torque/^ ^4TtR^G^^oj 

Ke)  ‘ 


is  twice  the  tensile  strain  e 
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measured  at  45°  with  respect  to  the  longitudinal  direction.  The 
proof  is  based  on  the  usual  Mohr's  circle  or  second  rank  tensor 
transformation  argument. 


Once  the  [A]  matrix  is  known,  the  composite  moduli  are  obtained 
from  the  system  of  four  independent  algebraic  equations  represented 
by  Eq  [C-3],  with  Eq  [C-7]  used.  The  four  equations  can  be  symme¬ 
trized  and  the  result  is 
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with  the  elements  of  the  symmetric  matrix  [M]  given  by 
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After  Eq  [C-13]  is  solved  for  the  composite  moduli  by  Grout  reduc¬ 
tion,  the  coupling  and  the  bending  stiffness  matrices  are  obtained 
from  Eq  [C-4]  and  [C-5] . 


C-13] 
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Two  comments  must  now  be  made.  First,  since  the  magnitude  of 
[B]  does  not  affect  the  determination  of  stiffness  and  moduli,  it 
seems  reasonable  to  assume  that  the  nonhomogeneity  of  [B]  will  not 
affect  the  analysis.  Second,  it  is  not  theoretically  necessary  to 
obtain  explicit  expressions  for  the  in -plane  stiffness  matrix  before 

c  c  c  c 

the  composite  moduli  are  computed.  The  moduli  ^22^  ^11’  ^12’  ^6 

may  be  obtained  directly  from  Eq  [C-9]  with  the  appropriate  loading 
conditions  specified  and  with  Eq  [C-1] ,  [C-3] ,  [C-4],  and  [C-7] 
used . 


Finally,  in  order  to  provide  a  comparison  with  theoretical  meth¬ 
ods  of  computing  moduli,  computations  have  been  made  for 
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Fig.  C-8  Photomicrograph  of  Weave  *’Line"  in  Cylinder 
1  (60X).  Surface  is  30  deg  with  Respect  to 
the  Axis  of  the  Cylinder.  Gage  10  is  also 
shown . 


NASA-Langley,  1965  CR— 266 


Fig.  C-6  Front  View  Just  Before 
Mandrel  has  Completed 
a  360-deg  Rotation 


Fig.  C-7  Rear  View  Just  Before 
Mandrel  has  Completed 
a  360-deg  Rotation 


"The  aeronautical  and  space  activities  of  the  United  States  shall  be 
conducted  so  as  to  contribute  ...  to  the  expansion  of  human  knowl¬ 
edge  of  phenomena  in  the  atmosphere  and  space.  The  Administration 
shall  provide  for  the  widest  practicable  ajid  appropriate  dissemination 
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